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Abstract 

The cosmological perturbations on a bouncing brane are studied. The brane is moving inside a 
Klebanov-Strassler throat where the infra-red region of the geometry is smoothly cut off. For an 
observer confined to the world- volume of the brane, this results in a non-singular bouncing mirage 
cosmology. We have calculated the scalar perturbations corresponding to the normal displacements 
of the brane. This is performed in the probe brane limit where the gravitational back-reaction of 
the brane on the bulk throat is absent. Our model provides a framework for studying the transfer 
of fluctuations from a contracting to an expanding phase. We find that the spectral index of the 
dominant mode of the metric fluctuation is un-changed, unlike what is obtained by gluing contract- 
ing to expanding Einstein universes with the help of the usual matching conditions. Assuming that 
the fluctuations start off in a vacuum state on sub-Hubble scales during the contracting phase, it 
is shown that the resulting spectral index on super-Hubble scales in the expanding phase has a 
large blue tilt. When the brane is moving slowly inside the throat and its kinetic energy is negligi- 
ble compared to its rest mass, one finds = 3. For a fast-rolling brane with a large kinetic energy, 
the spectral index is Ug = 2.3. This may put severe constraints on models of mirage cosmology. 
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I. INTRODUCTION 



Recent cosmological observations provide an increasingly clear picture of the current 
structure of the universe on large scales [T] . The universe is spatially flat to high accuracy, the 
background cosmology is homogeneous and isotropic, and there is a superimposed spectrum 
of almost scale- invariant, nearly Gaussian, and nearly adiabatic cosmological fluctuations. 
A viable theory of early universe cosmology must explain these observations. Inflationary 
cosmology [2] by far is the best model currently available which can explain these data. Fur- 
thermore, it successfully solves the flatness, homogeneity and horizon problems of Standard 
Big Bang (SBB) cosmology. On the other hand, current realizations of inflation are plagued 
by various conceptual problems [3] . Two of the problems which motivate this study are the 
presence of an initial singularity [4] which signals an incompleteness of the background cos- 
mology, and the "Trans-Planckian" problem, namely the fact that the fluctuations emerge 
from sub-Planckian wavelengths and hence from a zone of ignorance about the fundamental 
physics |31 E] • One should thus be open-minded towards alternative scenarios which may, 
in principle, be able to explain the data while avoiding the shortcomings associated with 
inflationary cosmology. 

Bouncing cosmologies can solve the horizon problem of the SBB cosmology. In bouncing 
cosmologies, each period of expansion is preceded by a period of contraction during which 
the co-moving Hubble radius is decreasing while the horizon continues to increase. Hence, 
in a bouncing universe the horizon can easily be made to be larger than the past light cone 
at the time of last scattering, the region over which the universe is seen to be homogeneous 
and isotropic. The entropy or size problem of SBB cosmology also disappears in the context 
of a bouncing cosmology: if the universe begins large in a phase of contraction, there is no 
reason why the initial entropy should be small. The flatness problem of SBB cosmology, 
however, persists (see e.g. [U]). 

A key question which determines the viability of any cosmological model with a bounce 
concerns the spectrum of cosmological perturbations. It is an extremely non-trivial challenge 
to flnd a scenario which is consistent with the current data. There are two issues involved. 
First, one must flnd a mechanism for producing an almost scale-invariant spectrum, and 
secondly, one must be able to evolve the spectrum through the bounce region. 

Recently, there have been several attempts to construct bouncing cosmologies, in particu- 
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lar the Pre-Big-Bang [7] and the Ekpyrotic [8] scenarios. Both, however, involve a singularity 
at the bounce point which prevents a consistent computation of the spectrum of fluctua- 
tions. (In a new version of the Ekpyrotic scenario, the singularity can be smoothed out 
by invoking ghost condensation [9J.) There have been attempts to construct non-singular 
bouncing cosmologies, making use of the interplay of curvature and matter with wrong sign 
kinetic terms [10], quintom matter pTI], k-essence [12] and higher derivative gravity actions 
[13] m]. An elegant way of obtaining a non-singular bouncing cosmology is by means of 
"mirage cosmology" [T3] . 

In this paper, we study the evolution of fluctuations through a non-singular bounce in 
the context of mirage cosmology. We focus mainly on two questions. First, does the spectral 
index of the cosmological fluctuations change when passing through this bounce? Secondly, 
does a scale-invariant spectrum emerge if we set up the fluctuations in a vacuum state early 
in the period of contraction. The answer to the second question is "no". Concerning the 
first question, we find that the spectrum goes through the bounce without change in the 
spectral index. This result is very interesting. In work on Pre-Big-Bang and Ekpyrotic 
cosmology where the contracting phase and the expanding phases were each described by 
dilaton or Einstein gravity and connected at the singularity by matching conditions [161 [IZj 
analogous to the Israel matching conditions [T8|, it was found [191 [20] that the spectral index 
changes: the dominant mode of the metric fluctuations in the contracting phase couples only 
to the decaying mode in the expanding phase. However, the applicability of these matching 
conditions has been challenged [21] and concrete studies [121 1221 122] have shown that it is 
possible that the spectral index does not change if the bounce is smooth (although other 
studies give differing results [131 [2l|). Mirage cosmology provides a simple setting to study 
the issue of the transfer of cosmological perturbations through a bounce in a clean way. 

II. THE BOUNCING BACKGROUND 

In "mirage cosmology" [T5| our universe is a probe D3-brane moving in extra dimensions. 
Due to the motion of the brane in the extra dimensions, the induced metric on the brane 
becomes time-dependent. For an observer confined to the brane this results in a cosmological 
expansion or contraction. In the probe brane limit, we neglect the back-reaction of the 
moving brane on the background. This is a good approximation when a single brane is 
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moving in the background of a large stack of branes or in the background created by a large 
number of flux quanta. 

The bouncing brane in a throat has been studied in [23], |2S], [21] and [2B]- The throat 
was taken to be the warped deformed conifold of the Klebanov-Strassler (KS) solution [29] , 
where the infra-read (IR) region of the geometry is smoothly cut off. The ultra-violet (UV) 
region of the throat is smoothly glued to the bulk of the Calabi-Yau (CY) manifold. The 
brane sets off from the UV region, moving towards the IR region. For an observer on the 
brane this corresponds to a contracting phase. At the tip of the throat, where the warp 
factor is stationary, the induced scale factor becomes stationary. This results in a bounce 
point. Finally, the brane bounces back to the UV region which corresponds to an expanding 
period. In order for the scenario to work, then one has to find a way to consistently connect 
this period to a late time Big Bang cosmology. 

One novel feature of this bouncing solution is that it is singularity-free. This is a mani- 
festation of the fact that the IR region is smoothly cut off. Furthermore, in and [22] it 
is shown that a cyclic cosmology can be achieved by turning on internal angular momenta. 
With angular momenta turned on, the brane motion is confined between IR and UV region 
of the throat which results in a cyclic solution. 

The metric of the background is given by 

ds^ =GMNdX^^dX^ =h{r)-^/^7]^^dxf'dx'' + g{r)dr^ + dsl. (1) 

Here, Gmn is the background metric, h{r) is the warp factor, r is the radial direction of the 
throat and ds^ is the metric along the internal angular directions. In what follows, we do 
not consider the angular motion of the brane and neglect the last term in Eq. ([T]). 

The space-time indices are represented by capital Latin letters M, N, ... while the indices 
used by the observer confined to the brane are represented by a, 6, .... The brane is equipped 
with coordinate In general, when the brane is moving non- uniformly, one cannot use the 
static gauge and {^"'} ^ {s^}. However, the fiuctuations along the spatial coordinates of the 
branes, are not physical and we can choose = for i = 1,2, 3. The time coordinate 
on the brane is defined by C,^ = t. In the homogeneous limit when the brane is moving 
uniformly and no perturbations are present on the brane world-volume, we can set X° = t. 

The action of the probe D3-brane is given by 

^ = 7^3 I rf'eV-|7afe|+/^3 / Am4))- (2) 
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Here, the first term is the usual Dirac-Born-Infeld (DBI) term and the second one is the 
Chern-Simons (CS) term. Also T3 and fi^ are the brane tension and charge, respectively, 
7ab is the induced metric on the brane, C(4) is the background Ramond-Ramond four-form 
field and P(C(4)) represents its pull-back onto the world volume Ai of the D3-brane. In the 
following, we consider the motion of a BPS brane and set /is = T3. 

The position of the brane in the extra dimension is parameterized by 

= X^iC) . (3) 
For a background space-time metric Gmn given in Eq. ([T]), we have 

lab - GmN-q^-q^ , /^(0(4)jil23 - ^O(4)o^23- W 

For the KS background [29], the radial part of the metric and the form of the warp factor 
is given by (for a detailed review of this background see [30] ) 

21/3 



where 



igsMa')Iiry/'K{r)-' , h{r) = {gMa'f2'/h-'/'l{r) , (5) 



./X /"°° , xcothx-1, , ,1/0 , (sinh 2r - 2r)^/^ , , 

/ r = / dx ^ sinh 2x - 2x^/3 , K(r) = ^ — — . 6 

^ ' Jr sinh^ X ^ J ^ K J 2V3 sinhr ^ ' 

In these expressions, Qs is the string coupling, M is the quantum number of the RR-fluxes 
turned on inside the at the bottom of the throat, and e'^^^ measures the size of this 
in the units of a'. The warp factor h{r) has a maximum at the bottom of throat r = 
and falls-off exponentially at large r. The function g{r) has a minimum at r = and 
increases like r^/^ for large r. The fact that h{r) and g{r) are non-singular and their first 
derivatives vanish at r = is crucial to obtain the bounce solution. In this background, the 
RR four-form field is given by 

C(4)oi23 = 97^h{ry^ ■ (7) 

At the homogeneous level the brane moves along the radial directions without any per- 
turbation, and, using Eqs. ^ and ([T]), we obtain 

S = J d'iC = T^ J d'ih-^ (-Vl-Zi^/^^r^ + l) , (8) 
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where the relation X*^ = t was used. Here and in what follows, an over-dot corresponds to 
a derivative with respect to t, the time measured by the observer on the brane. 

Integrating the radial Euler-Lagrange equation which follows from this action, one can 
construct a constant of integration 

E = fdC/dr-C. (9) 
This yields the following first order equation of motion for the position of brane r = R{t) 

R^ + V.MR) = , V.ffiR) = -^^^^l^^^^- (10) 

g{l + Eh) 

It is easy to see that E > and K// is always negative. Thus, the brane can move in 
the entire region < r < oo. On the other hand, to get a finite value for Newton constant 
(finite Mp), the KS throat should be glued smoothly to the bulk of the CY manifold at 
r = Tc- This implies that the brane would be moving in the region < r < Tg. We restrict 
our analysis to the throat region of the KS construction. 

An observer confined to the D3-brane feels cosmological expansion or contractions, de- 
pending on the direction of the brane's motion. For this observer, the four-dimensional 
metric is given by 

dsl = -dr"^ + a(r)2 rfx^ , (11) 

where x denote the three spatial coordinates and the proper time is defined by 

dr^ = -joo dt^ = h~^'^ (1 -gh^'^R^)de , (12) 

and the scale factor is given by a(r) = h {R{t)Y^^^ . 
The corresponding Friedmann equation is 

H^ = l-^^^[E^h^ + 2Eh)^ (13) 

where H = -^^^^^ is the Hubble constant and ' represents the derivative with respect to 
r. 

Fig. [l] shows how the scale factor behaves as a function of r and h as a function of r 
(in h{r) the prefactor has been set to 1). When the brane moves from the UV region to 
the IR region of the geometry, the observer sees a cosmic contraction. At r = 0, h' = 0, 
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FIG. 1: The warp factor h{r) and the scale factor a(r) are presented here. The scale factor is given 
by a(T) = h(R{t))^^^^. The contracting phase corresponds to the case when the brane moves from 
the UV region (large r) towards the IR region (small r). The bounce is at r = 0,t = 0, where 
H = 0. The expanding phase is the mirror image of the contracting phase. 



so H = 0. This corresponds to a bounce point Afterwards, the brane bounces back 
towards the UV region and the observer sees a cosmic expansion. The expansion rate H 
in the expanding phase is the mirror image of the expansion rate in the contracting phase. 
Using the asymptotic forms of the functions /(r) and K{r) given in Appendix B, one 
can see that immediately after the bounce, H enters a period of acceleration, followed by a 
period of deceleration. 

To perform the mode analysis in the following sections, it is very helpful to go to conformal 



time T), wheiea'^dr]'^ = dr"^. Using Eq. (10), we obtain 



III. FLUCTUATIONS OF THE BRANE POSITION 

In this section we study the fluctuation of the brane's position while it is moving inside 
the throat. We are particularly interested in the evolution of these perturbations measured 
by the observer confined to the brane. Our study is a generalization of what was previously 
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studied in |3T] . 

In general, the displacements of the brane can be decomposed into longitudinal and 
normal components. However, the longitudinal perturbations (which are along the tangential 
directions of the brane) are not physical and can be gauged away. Suppose is the space- 
like unit normal vector to the brane. Then the normal displacements of the brane are given 
by 

X^^ = X^^ + $n^. (15) 

Here, $ is the scalar field corresponding to the normal displacement and X^^ represents the 
brane position at the homogeneous background, i.e. the zeroth order brane position. It is 
important to note that just considering perturbation along the r direction is not consistent, 



at least when the probe brane is moving fast. As is clear from (19), the time-like component 
of the normal vector to the brane becomes important as the probe brane enters into the 
fast-moving regime (see Fig. [2]). 

At each point, we can span the 5-D space-time by the normal vector base {ea,n}, where 
{cq} are the tangent vector spanning the brane world volume. More, specifically 

ef = ^^X^ , 7a6 = G(e.,e,), (16) 

where G{ , ) represents the inner product defined by the background metric Q, and jab is 
the induced metric given in Eq. (0). The normal vector n^^ is defined by 



G{n,ea) = , G{n,n) = 1, (17) 

and the following projection relation holds 

r'X^a^^ = G^'"" - n^' . (18) 

For the background defined by Eq. ([T]), the non-zero components of the normal vector, 
and n'", are 



n 



° = g'/^h'/'R{l - h^'^gk^)-^/^ , = g-^'\l - h^/^gR^Y^'^ . (19) 



Using Eqs. (|4]) and ([T]), the action of the moving brane to all order in perturbations is 
S = -tJ d'^ (v^- ^q4)(X)] . (20) 
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FIG. 2: In this plot the hght cone for a fixed point on the brane is plotted. The dashed curved 
line represents the light-cone of the point when the brane is moving homogeneously governed by 



Eq. (10). The unit normal vectors are in — r plane and have components in both r and 



directions as is evident from Eq. (19). The curved line represents the light-cone when the normal 
perturbations Eq. (15) are present. Now n^'^ also has components in the tangential directions X*. 



The equation of motion for 6X^'^ which follows from this action is obtained in Appendix 
A, and is 

= (21) 

Here, K is the trace of the extrinsic curvature of the surface, K = 7"* Kab- The extrinsic 
curvature is defined by 

dX^^ dX^ 

= ef Dm = — — Dm , (22) 
where Dm is the five-dimensional covariant derivative which is compatible with the metric 



Gmn- One can check that at the homogeneous level Eq. (21 ) reproduces the solution studied 
in the previous section. 



To get the linearized equation of motion, we perturb Eq. (21 ) about the background ( 15 ). 
This is equivalent to expanding the action to second order in the perturbations. We relegate 
the details of the calculations leading to the perturbed equation of motion to Appendix A 
and here just give the final result 

- M"^^ = , (23) 



where is the Laplacian defined by the metric 'jat and the effective mass term is given by 

1 



n 



(24) 



In the above expression, Rmn is the background Ricci tensor. For the pure AdS background 
considered in [31] only the first two terms are present and the other two cancel out. 



We now calculate the different contributions to for our background. From Eq. (21) 
at the homogeneous level, we find 

h' 



K 



hg-y^c' = - 



(25) 



To calculate the other terms in M^, we need to express n^, n^,R and R as functions of r 



and E. To do this, we use Eq. (10) which allows R to be calculated. Plugging the value of 



R obtained this way into Eq. (19), one can show that 



n 



E^^h^/^ {2 + Ehf 



/2 



n 



g-'l' (1 + Eh) 



(26) 



Using these expressions for and ra'", then for the last term in we obtain 
a/-|7| R 



(27) 



For the term K°'^Kabi we have 



Agh^ 



DmUn D^' - Dm Dp n^'^ 



(28) 



To go from the second to the third line, the orthogonality condition in Eq. (17) and the 



projection relation Eq. (18) have been used. 



Finally, for the contribution from the Ricci tensor, we get 
-1 



Rmnu^u^ 



[3Eh{2 + Eh){2gh''^ - 2ghh' + hh'g') - 8ghh" + Ahh'g' + lOgh''^] (29) 



Plugging Eqs. (29), (28), (27) and (25) into the expression (24) for M we obtain 



E 
8g^ 



[h{2 + 3Eh) {g'h' - 2gh") + Agh'"^ ] 



(30) 



10 



IV. MODE ANALYSIS 



The differential equation (23) contains a Hubble friction term for $. As done in the 
usual theory of cosmological perturbations (see P2] for an in depth review and [33J for 
a pedagogical overview), we can extract the effects of the Hubble damping by the field 
redefinition 

0(r/,f) = a(r/)$(r7,f) = /i"^/^ $ , (31) 



where for convenience we use conformal time r]. Eq. (23) then becomes 

0-$:5^0+(/.-V^M^-^)0 = O. (32) 

i 

The above equation is analogous in structure to the equation which describes the evolution 
of cosmological perturbations in cosmology based on the usual four-dimensional Einstein 
gravity. On scales larger than the Hubble radius, the negative square mass term a^^/a is 
larger than the k"^ term coming from the spatial gradients. Thus, whereas on sub-Hubble 
scales the gradient term wins out and leads to the usual micro-physical oscillations, the 
oscillations freeze out approximately at Hubble radius crossing, and from then on the modes 
undergo squeezing. As is familiar from the usual theory of cosmological perturbations, scalar 
metric fluctuations are squeezed not with the factor a{ri) as would be the case if = 0, but 
with a modified factor which is usually called z{ri) which takes into account that ^ 0. 

It is important to point out that, as shown in [HT], for the observer on the brane our 
fluctuation $ acts as the usual Bardeen potential [M] which is generally denoted by \E'. 
In terms of in longitudinal gauge and in the absence of anisotropic stress, the metric 
including scalar metric fluctuations takes the form 

ds^ = -a^T]) [{1 + 2^)dr/2 _ (1 _ 2^)dx2] , (33) 
It is helpful to work in Fourier space where 

Mv) = [ d'x<P{r],x)e-'^'-^. (34) 



This transforms Eq. (35) into our desired mode equation 



+ {k' -Veff)(j)k = 0. (35) 



drj 
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FIG. 3: The plot on the left represents v^jf as a function of r. It has a global minimum at r = 
with Veff — —0.14 and a maximum at r ~ 1.9 with Veff — 0.04 and exponentially falls off for large 
r. The plot on the right represents Vgff as a function of p = (2Ef3)^^'^'q . The contracting and the 
expanding phases are mirror images of each other. 



where the "effective potential" is given by 

•Jeff 



a 



-h ' M i r 

Ah dri^ ^dri' \4: h 16 

' 2-^ + - , (36) 



8g \ h h g 



where in the final line Eq. (14) is used to replace d R/drj. 

We now make an observation which turns out to be important later on. A closer look 
at (36) reveals that in the vicinity of r = 0, Veff is always negative. This is easily inferred 



by noting that the 2nd and 3rd terms inside the bracket in (36) are both vanishingly small 
around r = to ensure the smoothness of the supergravity background. At the same time 
the fact that /i is a monotonically decreasing function as a function of r with a zero first 
derivative at r = guarantees that h" /h is negative. So one can conclude that around the 
bounce point, from the point of view of the brane observer, the effective potential seen by 
the brane fluctuation is negative. 

We try to obtain some analytical insights into the evolution of the mode functions given 



the effective potential of Eq. (36). Due to the complexity of the KS geometry involving the 
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functions /(r) and K{r), it is not possible to have an analytical solution for the entire range 
of r. However, in some limits we can perform analytical calculations which can capture 
qualitatively (and even quantitatively in the slow probe limit) the numerical results. 



Let us define h(r) = (31{r) and g = '-) I{rY^'^K{r) ^. Notice that Eq. (35) can be 
rewritten in the following form 



in terms of dimensionless variables 



+ -Veff)<pk = ^: (37) 



^^^^^ ' ^=^^^ ' ^e// = — ^t^e//. (38) 

The potential Veff is plotted in Fig. [3j As a function of the r-coordinate, it has a global 
minimum at the tip, r = 0, reaches a global maximum around r ~ 1.9, and falls off as 
rexp (— 2r) for large values of r. In the large r limit, using the asymptotic form of /(r) and 
K{r) given in Appendix B, one obtains 

^e// - ^^^^re--. (39) 

To solve the differential equation (37) analytically, we need to express Veff as a function 
of fj. This requires inverting r as a function of fj. The relation between fj and r is given in 



Eq. (14), which leads to 

r dr 

fj = ± . . (40) 

The negative branch corresponds to the case when the brane is falling towards the tip of the 
throat, fi = corresponds to the bounce point when the brane reaches the tip. The positive 
branch corresponds to the case when the brane bounces back towards the UV region of the 
throat. It is clear that Veff is symmetric under fj —fj. 

The above integral can not be computed for all value of r. However, in the limit of large 
values of r, when the KS solutions is well- approximated by pure AdS geometry, one can 
compute the integral analytically. Using the asymptotic forms of / and we have 

^2-1/63-1/2 / 



1 + 2Vi 4r'/3 



^2-1/63-1/2 r^/i --I - ^^'^ e^^/^) f4l1 
^ e ^l2'4'4' 3E/3 ^ > 
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where F{a, b; c, z) is the hypergeometric function. 



In two extreme regimes, Eq. (41) can be inverted. First consider the "slow-roll" limit 
when the brane is moving slowly and Ef3 « 1. This corresponds to the limit when the 
kinetic energy of the brane is much smaller than the brane rest mass in the DBI action in 
Eq. ([8|. Using the identity 



F(a, h- c, ^) = (1 - ^)-"F(a, c - 6; c, — -), (42) 

Z -L 

we obtain 

~ ± 3 X 2-^/6 e"/^ . (43) 



Using Eq. (39) for v^fj in the large r limit, the potential reads 

C In |r/| 

where C is an un-important coefficient (the reason will be clear later). 



(44) 



The other limit where Eq. (|41|) can be inverted is when E (3 e^'''^^ >> 1. This corresponds 



to the "fast-roll" limit, when the kinetic energy of the moving brane is comparable to the 



brane rest-mass. In this limit, using the identity given in Eq. (42) for hypergeometric 
functions, we obtain 

which in the large r limit leads to 
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It is interesting that in the fast-roll case, the f/-dependence of ^e//, up to a logarithmic 
correction, has the same form as in the inflationary cosmology. However, the constant 
prefactor is different: in the inflationary model the prefactor is 2, while for the fast-roll case, 
it is 10/9. The constant has its origin in the asymptotic limit of the KS geometry, where 
it reaches pure AdS space. As is well known in the context of inflationary or Ekpyrotic 
cosmology, the constant 2 is crucial to obtain a scale-invariant spectrum starting with the 
Bunch-Davis vacuum initial conditions. This change in the coefficient results in a non-scale- 
invariant prediction for the scalar spectral index n^, as we will see explicitly in the next 
section. 
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The effective potential v^fj as a function of p = {2E(3Y/'^fj is given in Fig. [sj Suppose 
tlie brane starts from the UV region heading towards the IR region of the throat, so r/ < 
initially. The modes start in the sub-Hubble region where —kfj » 1 and k"^ — Vejf > 0. 
They evolve as sub-horizon mode until "crossing the potential" (i.e. fc^ — Ve/f = 0) and 
then become super-Hubble. They remain super-Hubble until they cross the potential again 
near r = and become sub-Hubble. The evolution of the modes is sub-Hubble near the tip 
of the throat until the brane bounces back at r = or f/ = 0. The mode evolution for fj > 
is the mirror image of what happens for fj < 0. 

The scalar spectral index can be calculated at an arbitrary point like rf^, as long as the 
mode is super-Hubble but in the decelerating phase at late times. One may identify fj = rf^ 
with the surface of last scattering. 



V. SOLVING THE MODE EQUATION ANALYTICALLY 

Having obtained the effective potential Veff in the large r limit of both fast-roll and slow- 
roll cases, we can solve the mode equation and compute the spectrum of fluctuations using 
analytical approximations. First, recall that the scalar power spectral index Ug is given by 

pm = ^i0^r ~ k-^-'- (47) 

Let us start with the fast-roll case which is more amenable to analytic approximation and 
where the equations are more similar to those appearing in inflationary models. To solve 



the differential equation (37) analytically, we make a crude approximation here and neglect 



the logarithmic running in the effective potential (46). We simply replace the logarithm by 
one. Comparing with the exact result obtained numerically in the next section, we find that 
this is a reasonable approximation. 

We wish to compute the power spectrum of (pk at a late time 77* on super-Hubble scales. 
Let us denote the time that the mode k crosses the Hubble radius shortly after the bounce 
by ri2{k), the time it enters the Hubble radius shortly before the bounce by rii{k), and the 
time it exits the Hubble radius during the initial contracting phase by ?7o(fc). Due to the fact 
that the potential is very negative at the bounce point, in small k limit, the k-dependence 
of ri2{k) and rii{k) is negligible. It is crucial to note that the negativity of the potential in 
the vicinity of the bounce point is generic as it was argued in the previous sections. Also, 
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between rii{k) and ri2{k), i.e. in the bounce region, the mode is oscillating and its amplitude 
does not change. This implies that as long as the potential around the bounce point is 
negative, the bounce region does not introduce any k-dependence in the transfer matrix. 
Hence, the power spectrum of (pk at the late time 77=,, has the same spectral index as the 
power spectrum at pre-bounce time — r^*: 

P^{k,r],) = e\M^,)\^ = k'A\M-V*)\\ (48) 
where the amplification factor A is given by 

, ^ f \Mv*)\' \ f \Mv2{k))\' \ f IMvmi' \ 

\\Mv2{k)w) \\Mviik)WJ \\M-v*ikWJ ^ ^ 

and is, to a first approximation, independent of k. Thus, the spectral shape is determined 
by the growth of the mode functions between initial Hubble radius crossing at rio{k) and the 
time — 77*: 

P,ik,v.) - fc^ '|^^|"Jj|jV ^(^o)r (50) 
In the fast-roll case, the dominant solution of (pk in the contracting phase is given by 

0fc(^)~r'/'- (51) 
Since rjo^k) is given by fc^ = Veff{rio{k)) we have 

r]o{k) ^ k~K (52) 



Inserting (51) and (52) into (50) and assuming that the modes start out on sub-Hubble 
scales during the period of contraction (i.e. for rj < r]o{k)) in the vacuum state with the 
Bunch-Davis initial condition 

e-''^\ (53) 



/2k 

we obtain 



P^k) (54) 

and thus the spectral index is predicted to be = 5/3. 

To be slightly more precise, using the large r limit of Veff form Eq. (46), the solution of 



(|37|) is given in terms of Hankel functions of first and second kind 

A,Hy,'-'\-kfl)+A2Hy,^^\-~kfj)] , (55) 
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where Ai and A2 are two constants of integration. Like in inflationary models, only the 
Hankel function of first type matches with the initial vacuum condition since 



Hj/e^'\x » 1] 




TTX 




(56) 



which implies that A2 = 0. Thus, we obtain the scaling of 0^ with 77 given in (51), and, as 
shown above, this leads to 



Hg = 5/3 . 



(57) 



Interestingly enough, the spectrum is not scale invariant. As explained before, this is due 



to the numerical factor 10/9 in Eq. (46) for "We//- This particular number originated from 
the fact that the UV region of the throat is well-approximated by pure AdS geometry. This 
indicates that for pure AdS geometry, the spectral index will be = 5/3. As explained 
before, we have neglected the logarithmic correction into "Ce//, which originates from the 
logarithmic correction to the AdS geometry in the KS solution. Considering the logarithmic 
correction, our numerical analysis (presented in next section) show that Ug ~ 2.3. 

The case of slow-roll proved to be more involved. Even when we neglect the logarithmic 



correction for Ve/j in Eq. (44), we could not solve the differential equation (37) simultane- 
ously for both sub- Hubble and super- Hubble modes. However, we can solve it for sub-Hubble 
and super-Hubble modes separately. For sub-Hubble modes the solution is as given by vac- 



uum state Eq. (53). For super-Hubble modes , the solution of Eq. (37) with effective 



potential Eq. (44) is given in terms of Bessel and Neumann functions 

5,J_v4(^) + 52F-i/4(^; 



(58) 



At the transition fjQ{k) point when k'^ — Ve/f = 0, this solution should smoothly go over to 
the solution for the sub-Hubble modes. The transition point is given by 



(59) 



Since we are interested in large scale perturbations, k « 1, so ?7o(/c) >> 1 and we can use 
the small argument approximation for the Bessel and Neumann functions 



Jy{x « 1) 



X 



Y,{x « 1) 



-X 



(60) 
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Using these approximations in Eq. (58 



), we obtain 



Mm) ^ (fii + fi2^o), (61) 

where Bi and B2 are new constant of integrations related to Bi and -82- Matching this 



solution to the sub-Hubble solution Eq. (53) at r]Q(k), we find 



Bi ~ -^(l + ikfio(k))e''^^°^''^ , B2 ~ ^^Vke-'''^°^''K (62) 
For large scale perturbations where k ^ 0, Bi scales like while B2 scales like k^^"^ and 



from Eq (61) we obtain (pkivo) ~ ^ ^^'^ ■ Using this in Eq (50), one obtains 



Hs = 3 . (63) 

We have calculated the spectral index numerically and found = 3 ± 10^^, in agreement 
with the above analytical value. 

Indeed, one can argue that for potentials of the form Veff ~ with n > 2, the spectral 
index is always equal to 3. On the other hand, n = 2 is the critical value where the spectral 
index crucially depends on the numerical coefficient in Se//- For a potential of the form 
Veff = cr/~^, one can show that = 4 — Vl + 4c. For infiationary models in the limit 
where slow-roll corrections are ignored c = 2 and one obtains = 1. In our fast-roll case 
c = 10/9, which leads to Ug = |, as we have shown before. 

VI. NUMERICS 

In this section we discuss our numerical results. In order to solve the equation of motion 



for (j)k (37) numerically, we find it more convenient to work with the coordinate r related 



to f] through (40). In this coordinate, we look for a left moving wave solution initialized at 



very large values of r which eventually refiects back to the large r region. The equation of 



motion is initialized with the Bunch-Davis vacuum initial data given in Eq (53). Boundary 
conditions at the origin of the r-coordinate are easily found to be the impact boundary 
conditions 



(Pkin\r=0+ (Pkout\r=0~j 

''4>kin, _ d(j)kout, 

—. |r=0+ — ~, |r=0- 

dr dr 
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(64) 




-6.5 -6.0 -5.5 -9.30 -9.28 -9.26 -9.24 



log(k) log(k) 

FIG. 4: The /c-dependence of (pkiv*) is numerically plotted. The figure on the left corresponds 
to the slow-roll limit with the slope equal to —0.5 and Hs = 3, while that of the right hand side 
corresponds to the fast-roll case with slope equal to —0.85 and Ug = 2.3. The data points are 
shown by red triangles. 



where (pkin and (pkout refer to the incoming and outgoing waves at the origin, respectively. 



The quantity of interest to us is the spectral index defined by (47). As usual, (pk is evaluated 
at a convenient point in the super-Hubble region. It is worth emphasizing that picking a 
different point will only amount to a change in the amplitude of the power spectrum and 
will not alter its k-dependence or the index as long as it lies within the super-Hubble region. 



We solve (37) for different values of k and read off the index from the log-log plot of A; |0fc| 
versus k, generated numerically for a range of k values. We observe a negligible running of 
the index, for either small or large values of the parameter EjS, as we dial k. We made sure 
that a change in the location where the initial Bunch-Davis vacuum condition was imposed 
on our numerical solution did not influence the resulting index. Fig. |4] illustrates our results 
for the index in two different limits for E(3. The left side represents the results in the small 
Ef3 limit. The right side gives the result in the opposite limit where EjS is large (by which we 



simply mean the limit where the second term under the square root in (40) can be ignored). 
Note that the limit itself depends on the wavelength of the observed large scale structures, 
k~f^g, around which the index is calculated. 
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FIG. 5: The wave function is presented here. The red(blue) curve corresponds to contract- 
ing(expanding) phase. As time increase, i.e. fj increases, the wave function also increases. This 
indicates that the dominant mode in the contracting phase is matched to the dominant mode in 
the expanding phase. 



Fig. [sj shows a sample wavefunction (its modulus to be exact) plotted for k = 0.0034. 
The lower curve represents the wavefunction in the contracting phase. The amplitude is 
growing as the bounce point r = is approached. The upper curve is the amplitude of the 
wavefunction in the expanding phase. As is apparent, the amplitude is growing in time (i.e. 
as fi increases). This demonstrates explicitly that the dominant mode of (pk in the contracting 
phase couples with un-suppressed amplitude to the dominant mode in the expanding phase, 
unlike what is obtained in models in which a contracting Einstein cosmology is matched 
to an expanding Einstein cosmology at a singular hypersurface making use of the usual 
matching conditions of [161 E] • 

Note that in the simulation of Fig. 4, the initial condition was imposed at r = 10, where 
the brane is deep in the UV region, well-approximated by the AdS geometry. 

VII. BACKREACTION 

As discussed in [23], we need to check that the gravitational instability associated with 
long wavelength gravitational perturbations does not destabilize the homogeneous back- 



20 



ground. The length scale associated with this instability is the Jeans length given by 
Lj = Vgj y/G~p where Vg is the speed of sound, G is the Newton's constant and p is the 
energy density on the brane. The time scale associated with this instability is tins > Lj. 
For our model, from Eq. (jsj), we obtain 



p = T^h-^ (1 - x/l - gh^H^) 
1 + Eh' 

In the slow-roll limit, one can see that p ~ ET3 is constant. 
Using G < gill for our compactification [25], we have 



(65) 



tins > Lj = -j== > Vshy — ^ — • (66) 

Suppose one releases the brane from the point r = r^, with the initial velocity v^. Since 
the brane is accelerating towards the tip of the throat, the time for the bounce satisfies 

^bounce — ; (^') 

where is the physical distance traveled by the brane before reaching the bottom of the 
throat. One can show that 

V, = ^/Eh{n){2 + Eh{r,)) (68) 

Discarding factors of order unity, taking f ^ ~ 1 and using the fact that g4/3g2r,/3 ^ ^2 
can show 

, 3/4 

^ < -j^[{l + Eh{nm + Ehin))]-'/\ (69) 

In the slow-roll limit when Eh{r^) « 1, one reaches the bound obtained in [23]. We see 
that this bound gets even weaker for the fast-roll limit, when Eh{r^) » 1. As an estimate, 
one can take r,, ~ (5 — 10), when the brane is deep inside the throat. Also, to trust the 
supergravity limit, we take gsM » 1. In this limit it is thus justified to neglect the effect 
of gravitational instability on the brane. 

VIII. DISCUSSION 

We have studied the evolution of cosmological fluctuations in a mirage cosmology setup 
in which the observer lives on a BPS D3-brane which is moving into and out of a Klebanov- 
Strassler throat. This provides a simple model of a non-singular bouncing cosmology. 
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Our main results are twofold. First, we find that the growing mode of the cosmological 
fluctuations in the contracting phase couples without suppression to the growing mode in the 
expanding phase. This is unlike what happens in a setup in which the fluctuations are treated 
in Einstein gravity in both the expanding and contracting phase, and then matched through 
a singular hypersurface using the analog of the Israel matching conditions. Negativity of 
the effective potential felt by the perturbations in the neighborhood of the bounce point 
appears to be generic in the models of the kind studied in this paper. As we saw, this leads 
to the fact that the bounce region does not contribute any ^-dependence to the spectral 
index. Secondly, we find that if we set off the modes in their Bunch-Davies vacuum on sub- 
Hubble scales in the contracting phase, a final spectrum in the expanding phase which is not 
consistent with observations will emerge. The specific spectral index depends on whether 
the brane is moving fast or slow. In the latter case, the resulting spectral index is Us = 3, in 
the former is closer to 2 but still completely inconsistent with the data which demands 
Us = 0.95 ±0.05. 

In this analysis we did not take into account the effects of volume modulus stabilization, 
such as can be achieved by wrapped D7-branes, on the mobile D3-brane. In a realistic 
model where all back-reactions are included, the brane feel an attractive force towards the 
tip of the throat [36l [37] . This changes Vef / and it is interesting to see whether or not these 
effects can modify the spectral index. It is also possible that mirage cosmology motion in a 
different type of throat might lead to a spectrum consistent with the cosmological data. It 
is also possible that a pre-cursor phase which involves extra physics (as assumed e.g. in the 
Ekpyrotic scenario) leads to a consistent spectrum. Finally, it is possible that, if the bounce 
phase lasts sufficiently long, thermal fluctuations of a gas of closed strings with winding 
modes will generate a scale-invariant spectrum [38l |39| HQ]. 
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APPENDIX A: LINEAR EQUATION OF MOTION 



Here we present the details of the calculations leading to the linear perturbation equation 



(23). Although we perform the analysis for a D3-brane moving in a five dimensional back- 



ground but the result is applicable to a Dp-brane with co-dimension one, i.e. p = D — 1, 
where D is the background space-time dimension. 

As was briefly described in Section 2, at each point in space-time we can choose the 
normal vector base {ea,n}, such that 



ef = , lab = G(e,,e,) 



and 



G{n,ea) = , G{n,n) 

where G{u,v) = Gmnu^v^ ■ 

The following projection relation also holds: 



(Al) 



(A2) 



MN 



We are interested in the normal displacements of the brane 



M 



(A3) 



(A4) 



Under this transformation, we have 



Sjab = ^{GMN,pn''X^X^ + 2GMNn'^X^) 



(A5) 



where Kab is the extrinsic curvature of the surface defined in Eq. (22 



The equation of motion for the normal displacement $ from the action (20) is obtained 



by 



= -T, j <Fi 



95A° 



(A6) 
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which produces Eq. (21). 



To get the hnearized equation of motion we perturb Eq. (21) around Eq. (A4). The 



perturbed equation of motion is to some extent involved. Here we present the outhne of the 
calculations. The linear equation is obtained by perturbing the left hand side (LHS) and 
right hand side (RHS) of Eq. (21) separately. 
For the LHS, we have 

6iLHS) = 5( v^^) = [(K' - 2K'''Kab)^ + ^'SKab] ■ (A7) 



To obtain the terms proportional to $, the relation (A5) for S'jab was used. To calculate 



6Kab, it is very useful to perform the analysis in locally Gaussian normal coordinates (GNC), 



where at each point G 



MN.P 



QS 



0, where Tq^ is the connection compatible with the 



metric Gmn- Starting with the expression for Kah given in Eq. (A5), in GNC we have 



^'^'dK,, = -n^rfi (G^^^ - n*^n^)GM^,PQ + I'^'G mn{5X^! n'l + Xpn';'^) (A8) 



where to obtain the terms proportional to $ the projection identity Eq. (A3) was used. We 



now calculate the last two terms in Eq. (A8) separately. To do that we note that 



<5n^^^ = A^e'^ +Bn 



M 



(A9) 



where 



B 



^ MNP 



^ -^-^^,-^-^^GMN,pX'ln''^n^ + GMNn^'n^b)^ 



Using this for the third term in Eq. (A8) we obtain 



7 '^MNO^^b 



r''GMNn'':n'l^ 



ab ■ 



(AlO) 



(All) 



To go from the first line to the second line, the relation n^^nM,a = which holds in GNC 
was used. The final identity in the above equation is easy to prove in GNC following a chain 
identity for partial derivatives. 
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Similarly for the last term in Eq. (A8) we have 



= jT^ + A^^'CmnX^X^. (A12) 

One notes that the term containing B vanishes in GNC while the term containing B^a 
vanishes due to relation G{n, Ca) = 0. 

After some algebra, and working in GNC, one can show that 

- l'''^,ac - l'::^,a (A13) 

and 

A'^Y'GmnX^, X,f, = -r'l^'GMNX^, . ( A14) 

Combining all the terms, we obtain 

7'^'Gm7vX^5< = -V^^-K'^'Ka,^-i'''GMNn'''nl^^ 

- n^n^) n^'n'^GMN,PQ $ (A15) 

where is the Laplacian defined by the metric 7afe. 



Combining Eqs. (A15), (All) and (A8) we obtain 



Wn'^n'-n^gPQ^MN^ - ^'GuNn'^'nl,^ . (A16) 

We need to get rid of the last term above containing second derivatives of n. To do this, 
we take derivatives of the identity G(n, n) = 1 in GNC, which leads to 

^'^"GMNU^'nl, = -^n^n^ (G^^ - n'^)G m n ,pq - K'^'K^b • (A17) 



Using the above expression in Eq. (A16), we have 

j'^'SKab = -V'<l> + K'^'Kab - RMNn^n'' , (A18) 
where Rmn is the Ricci tensor of the background geometry. 
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Thus, our final expression for S{LHS) in Eq. (A7) is 
6{LHS) = 



(A19) 



in agreement witli 



For tlie RHS of Eq. (21) we have 



6{RHS) = C['4)n"K-i?n°) + q4)KnO-nV) ^ + C[^){6n'' - R 5n^) ^ . 



(A20) 



Using the formula (A9) for Sn, one can show that 

(Goo + GrrR'^) 



6n' -R6n^ 



2Goo drr {n^R — 



n 



R 



(A21) 



where to get the final line the relations G{n, n) = 1 and G{n, e^) = were used. 

tain 



Combining Eqs. (A21) and (A20) we obtain 



6{RHS) 



R 



(A22) 



Combining Eqs. (A22) and (AT9|, the desired linear equation of motion from Eq. 

1 



(21) is 



^("»-|)(qV'0,, 



$ = 0. (A23) 



APPENDIX B: APPROXIMATIONS FOR /(r) AND K{r) 

The following approximation formulae for the KS background are useful: 

K{r^O) (2/3)^/3 + C(r2) , 
K(r^oo) ^ 2i/3g-^r/3^ 

/(r^O) ^ 0.72 + O(r2), 
I(r ^oo) ^ 3.2-i/Mr-^')e-^'^/3_ 



fBll 
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